An efficient meshfree point collocation moving least squares method to solve the interface problems with nonhomogeneous jump conditions

We are going to study a simple and effective method for the numerical solution of the closed interface boundary value problem with both discontinuities in the solution and its derivatives. It uses a strong‐form meshfree method based on the moving least squares (MLS) approximation. In this method, for the solution of elliptic equation, the second‐order derivatives of the shape functions are needed in constructing the global stiffness matrix. It is well‐known that the calculation of full derivatives of the MLS approximation, especially in high dimensions, is quite costly. In the current work, we apply the diffuse derivatives using an efficient technique. In this technique, we calculate the higher‐order derivatives using the approximation of lower‐order derivatives, instead of calculating directly derivatives. This technique can improve the accuracy of meshfree point collocation method for interface problems with nonhomogeneous jump conditions and can efficiently estimate diffuse derivatives of second‐ and higher‐orders using only linear basis functions. To introduce the appropriate discontinuous shape functions in the vicinity of interface, we choose the visibility criterion method that modifies the support of weight function in MLS approximation and leads to an efficient computational procedure for the solution of closed interface problems. The proposed method is applied for elliptic and biharmonic interface problems. For the biharmonic equation, we use a mixed scheme, which replaces this equation by a coupled elliptic system. Also the application of the present method to elasticity equation with discontinuities in the coefficients across a closed interface has been provided. Representative numerical examples demonstrate the accuracy and robustness of the proposed methodology for the closed interface problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1031–1053, 2015     

Moving boundary truncated grid method: Multidimensional quantum dynamics

The moving boundary truncated grid (TG) method is used to study wave packet dynamics of multidimensional quantum systems. As time evolves, appropriate Eulerian grid points required for propagating a wave packet are activated and deactivated with no advance information about the dynamics. This method is applied to the Henon-Heiles potential and wave packet barrier scattering in two, three, and four dimensions. Computational results demonstrate that the TG method not only leads to a great reduction in the number of grid points needed to perform accurate calculations but also is computationally more efficient than the full grid calculations.     

大气环境下的钢筋混凝土桥梁时变可靠度分析

为了探究钢筋混凝土桥梁抗力退化对桥梁可靠性能的影响,综合考虑一般大气环境下混凝土强度、钢筋强度、混凝土与钢筋黏结性能的降低,建立普通钢筋混凝土桥梁的抗力退化时变模型,分析规范车辆荷载作用的桥梁时变可靠度。研究结果表明,(1)一般大气环境下,桥梁抗力随着服役时间的增加而减小,其中钢筋强度退化对桥梁抗力的影响最大,其次是钢筋与混凝土黏结性能降低,而混凝土强度退化对桥梁抗力的影响相对较小; (2)桥梁可靠指标随着服役时间的增加而减小,密集行车占比越大,桥梁承载能力失效的概率就越大;非平稳车载作用下的桥梁可靠指标小于平稳车载,桥梁需要维修补强的时间比平稳车载提前; (3)多梁式简支梁桥在建造时应该提高边梁的安全储备,在管养时可以采取限载措施进行干预。     

阻尼对空爆荷载等效静载动力系数的影响

为考查阻尼参数对空爆荷载等效静载动力系数的影响,理论推导了空爆荷载下结构等效单自由体系弹塑性位移解及延性比解,设计并计算了阻尼比0.000 1～0.1、延性比1～4的20种典型工况的动力系数,并与现行抗爆设计规范动力系数公式结果进行了对比。结果表明:阻尼比小于0.000 1时可基本代表无阻尼状态,阻尼比0.01的动力系数比无阻尼的最大降低幅度为2.08%,数值差异很小,因此阻尼比为0.01以内时,可忽略阻尼对动力系数的影响;阻尼比0.05的动力系数比无阻尼的降低幅度约9.92%,数值差异较大,认为阻尼比0.05以上时将具有明显的经济效益;现行设计规范动力系数更适用于柔性结构体系,运用于刚性结构抗爆设计时,计算误差较大,对阻尼比较小的结构设计更不利。     

分段滤波用于消除脉搏血氧检测中多种运动干扰的研究

脉搏血氧饱和度检测仪可以无损、实时地检测动脉血液中的含氧百分比,从而快速对人体的呼吸状况以及心肺功能作出判断,但其检测精度极易受到测试部位运动的影响.人体测试部位的运动大致可以分为突发性运动干扰以及周期性的运动干扰.突发性干扰表现在局部脉搏波的波形突变上,而周期性的干扰体现在多个局部脉搏波的波形的周期性变化上.文章有针对性地提出了一种将微分阈值分段与迭带中值滤波相结合的方法来剔除突发扰动的方法,并且在上诉方法上再加上一级多点移动平滑以消除周期性运动干扰.实验证明该方法在消除多种运动伪差上有着较好的效果,且该方法计算简单,易于实现,有着一定的实际应用价值.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

Modified characteristic finite difference fractional step method for moving boundary value problem of percolation coupled system

For the coupled system with moving boundary values of multilayer dynamicsof fluids in porous media,a characteristic finite difference fractional step scheme appli-cable to the parallel arithmetic is put forward.Some techniques,such as the change ofregions,the calculus of variations,the piecewise threefold quadratic interpolation,themultiplicative commutation rule of difference operators,the decomposition of high orderdifference operators,and the prior estimates,are adopted.The optimal order estimatesin the l2norm are derived to determine the error in the approximate solution.This nu-merical method has been successfully used to simulate the flow of migration-accumulationof the multilayer percolation coupled system.Some numerical results are well illustratedin this paper.     

Development of a multi‐scale ocean model by using particle Laplacian method for anisotropic mass transfer

The direct injection of CO₂ into the deep ocean is one of the feasible ways for the mitigation of the global warming, although there is a concern about its environmental impact near the injection point. To minimize its biological impact, it is necessary to make CO₂ disperse as quickly as possible, and it is said that injection with a pipe towed by a moving ship is effective for this purpose. Because the injection ship moves over a spatial scale of O(10²km), a mesoscale model is necessary to analyse the dispersion of CO₂. At the same time, since it is important to investigate high CO₂ concentration near the injection point, a small‐scale model is also required. Therefore, in this study, a numerical model was developed to analyse CO₂ dispersion in the deep ocean by using a fixed mesoscale and a moving small‐scale grid systems, the latter of which is nested and moves in the former along the trajectory of the moving ship. To overcome the artificial diffusion of mass concentration at the interface of the two different grid systems and to keep its spatial accuracy almost the same as that in the small‐scale, a particle Laplacian method was adopted and newly modified for anisotropic diffusion in the ocean. Copyright © 2010 John Wiley & Sons, Ltd.     

A study of gridless method with dynamic clouds of points for solving unsteady CFD problems in aerodynamics

When solving unsteady computational fluid dynamics problems in aerodynamics with a gridless method, a cloud of points is usually required to be regenerated due to its accommodation to moving boundaries. In order to handle this problem conveniently, a fast dynamic cloud method based on Delaunay graph mapping strategy is proposed in this paper. A dynamic cloud method makes use of algebraic mapping principles and therefore points can be accurately redistributed in the flow field without any iteration. In this way, the structure of the gridless clouds is not necessarily changed so that the clouds regeneration can be avoided successfully. The spatial derivatives of the mathematical modeling of the flow are directly determined by using weighted least‐squares method in each cloud of points, and then numerical fluxes can be obtained. A dual time‐stepping method is further implemented to advance the two‐dimensional Euler equations in arbitrary Lagarangian–Eulerian formulation in time. Finally, unsteady transonic flows over two different oscillating airfoils are simulated with the above method and results obtained are in good agreement with the experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

A ghost fluid Lattice Boltzmann method for complex geometries

A ghost fluid Lattice Boltzmann method (GF‐LBM) is developed in this study to represent complex boundaries in Lattice Boltzmann simulations of fluid flows. Velocity and density values at the ghost points are extrapolated from the fluid interior and domain boundary via obtaining image points along the boundary normal inside the fluid domain. A general bilinear interpolation algorithm is used to obtain values at image points which are then extrapolated to ghost nodes thus satisfying hydrodynamic boundary conditions. The method ensures no‐penetration and no‐slip conditions at the boundaries. Equilibrium distribution functions at the ghost points are computed using the extrapolated values of the hydrodynamic variables, while non‐equilibrium distribution functions are extrapolated from the interior nodes. The method developed is general, and is capable of prescribing Dirichlet as well as Neumann boundary conditions for pressure and velocity. Consistency and second‐order accuracy of the method are established by running three test problems including cylindrical Couette flow, flow between eccentric rotating cylinders and flow over a cylinder in a confined channel. Copyright © 2011 John Wiley & Sons, Ltd.